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We report on a study of the exchange tunneling reaction D+HD—:-D 2 +H in a pure solid HD 
matrix and in a D 2 matrix with 0.23% HD admixture at temperatures between 130 mK and 1.5 K. 

We found that the exchange reaction rates, k^xHO ~ 3 x 10~^^cm®s“^ in the pure HD matrix, and 
kexD 2 = 9(4) X 10“^®cm“® in the D 2 matrix are nearly independent of temperature within this 
range. This confirms quantum tunnelling nature of these reactions, and their ability to proceed 
at temperatures down to absolute zero. Based on these observations we concluded that exchange 
tunneling reaction H+H 2 —>'H 2 +H should also proceed in a H 2 matrix at the lowest temperatures. 

On contrary, the recombination of H atoms in solid H 2 and D atoms in solid D 2 is substantially 
suppressed at the lowest temperatures as a result of increasing of violation for resonance tunneling 
of atoms when they approach each other. 


I. INTRODUCTION 

Solid hydrogen and deuterium form a special class of 
quantum solids, where due to a small mass and weak 
intermolecular interactions, the effects of quantum tun¬ 
neling play an important role. Light atomic impurities 
like H and D stabilized in matrices of the solid hydrogens 
are able to migrate from one lattice site to another by a 
repetition of tunneling reactions^: 

H-bH2^H2-bH. (1) 

D + D 2 -T- D 2 + D (2) 

A clear consequence of this is that a recombination of H 
and D occurs when two atoms encounter each other in 
neighboring lattice sites. Tunneling reactions also take 
place in H 2 — D 2 mixtures and involve both hydrogen 
isotopes: 

D -b H 2 ^ HD -b H. (3) 

D -b HD ^ D 2 -b H (4) 

A spectacular production of H atoms occurs as D atoms 
combine with H 2 in reaction ([31) and HD in reaction ([3|). 
Reactions (|4|) and (|3|) are well-known in the gas phase at 
higher termeratures and were intensively studied in H 2 - 
D 2 solids |2|-|^ at temperatures down to ^1 K. The differ¬ 
ence in the zero-point energies between the products and 
reactants (about 400 K) results in a preferential creation 
of H atoms while the reverse reactions are endothermic 
and do not proceed at low enough temperatures. The 
rates for the tunneling reactions (| 3 |) and were theo¬ 
retically calculated by Takayanagi et al. Q and Hancock 


et al. The rate of the reaction 0 was measured 

experimentally in the temperature range 1.9-6.5 K in a 
series of works by Miyazaki and Kumada 0 , 0 . Reaction 
m proceeds two orders of magnitude faster audits rate 
was measured only recently by Kumada 0|. 

The main signatures of a tunneling reaction are weak 
dependence of its rate on temperature and a large iso¬ 
topic effect. For a single atomic impurity in a perfect 
crystal tunnelling motion due to the reactions © and 
[ 2 ] may be very fast because of resonance in the energies 
of the initial and final states, since all locations for the 
atoms in the lattice are equivalent to each other. Kagan 
et al. [13 pointed out that crystal defects and irregu¬ 
larities violate the condition of resonant tunneling of H 
atoms in H 2 matrix which affects the rate of a tunneling 
reaction. Kumada et al. showed that the ortho¬ 

molecules as 2nd and 3rd nearest neighbors significantly 
decrease the rate of the atomic hydrogen tunneling reac¬ 
tion in a para-H 2 matrix. According to Kagan, tunneling 
in irregular solids may be stimulated by phonons, which 
force tunneling rates to be proportional to T for a single¬ 
phonon direct process or to T"=5-9 for a two-phonon 
Raman process 03 . It was found experimentally that 
the recombination rate of H atoms in solid H 2 follows 
a linear dependence on temperature in the temperature 
range 1.2-4K [^, while the rate of the exchange reaction 
eq. ([4|) remained nea^ temperature independent within 
experimental errors |8|. Significant divergence from the 
linear law was observed for the H recombination rates 
in solid H 2 upon reducing temperature below 1 K where 
the H atom recombination decreased by more than two 
orders of magnitude for samples cooled down to 150 mK 

MM- 

Since in solid D 2 with a small H 2 or HD admixture atomic 
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impurities need to approach a H 2 or HD molecule before 
taking part in reactions ([3]) and o, these reactions as 
well as the recombination of atoms are two stage pro¬ 
cesses: 1) diffusion due to tumrelliirg exchange followed 
by 2) the reaction itself. The matrix of pure HD repre¬ 
sents a different case, since the atomic impurities in the 
lattice are always surrounded by 12 HD neighbours. The 
diffusioir process is not needed for the reaction ([4]) in a 
solid HD matrix. 

In the present work we report on an experimental study 
of the isotopic exchange reaction corresponding to eq. 
in a HD matrix in a temperature range 0.13-1.5 K. An 
opposite limit where D atoms need to diffuse in order to 
encounter a HD molecule to react with was studied in a 
solid mixture of D2:0.23% HD. Our cryogenic system is 
based on a dilution refrigerator, which allows continuous 
measurements of the reaction kinetics for experiment du¬ 
rations of the order of weeks and months, substantially 
longer than in any previous work. This makes possible 
measurements of very slow chemical processes with very 
small concentrations of reagents, not accessible before. 
Our data for the rate of the reaction Q at T=0.13-1.5 K 
in a pure HD matrix agree well with previous works [J @| 
and confirm a very weak temperature dependence in this 
range. We found that the reaction ([4]) in the D 2 matrix is 
also temperature independent, which indicates that the 
non-resonant diffusion due to the exchange reaction ([5]) 
is not a limiting stage. The rate of the reaction (|3]) in the 
D 2 matrix proceeds about 3 times slower than in a pure 
HD matrix. At temperature of 150 mK we have not ob¬ 
served any recombination of D atoms in the D 2 matrix, 
while the exchange reaction o was well visible. Since 
the recombination also involves the diffusive process, we 
anticipate that the atomic motion due to the reaction ([J) 
is slowed down substantially when the D atoms approach 
each other. Based on that hypothesis we concluded that 
in the previous work the exchange tunneling chemical re¬ 
action H+H 2 —>'H 2 +H also takes place in solid molecular 
hydrogen at 150 mK while the atomic recombination is 
suppressed by non-resonant tunnelling when two H atoms 
apporoach each other (T^. 

II. EXPERIMENTAL SETUP 

Experiments were carried out in the sample cell (SC) 
shown in Fig. mEl. The SC is located in the center of a 
4.6 T superconducting magnet and is cooled by a dilution 
refrigerator down to about 100 mK. The main diagnostic 
tools in our experiments are an electron spin resonance 
(ESR) spectrometer employing a Fabry-Perot resonator 
operating at 128 GHz frequency 0 and a quartz mi¬ 
crobalance (QM) to measure the sample thickness. ESR 
enables measurements of H and D concentrations and 
hyperfine state populations. Absolute atomic concentra¬ 
tion calibration is obtained from the dipolar broadening 


to ESR 



Figure 1. Schematic of the sample cell. QM: Quartz crys¬ 
tal microbalance (the bottom mirror of the Fabry-Perot res¬ 
onator), HNMR: 910 MHz helical resonator for RF discharge 
and H NMR. 

of the ESR spectra 0 or calorimetrically by recording 
gas phase hydrogen recombination via measuring heat re¬ 
leased due to reactions 00. Solid films of hydrogen 
isotopes were deposited on the flat ESR resonator mirror 
which also serves as the QM electrode. 

The samples were condensed directly from a room tem¬ 
perature reservoir through an electrically heated capil¬ 
lary (Fig[T]). The capillary end is directed at the QM 
which is kept at temperature of 0.5-1 K. The QM makes 
it possible to measure the thickness of the growing films 
starting from 0.2 monolayer and to use small deposition 
rates to decrease overall heating due to the sample prepa¬ 
ration. After finishing preparation of a solid film on the 
surface of the QM, H and D atoms were created by dis¬ 
sociation of HD or D 2 molecules in situ in the solid by 
electrons created by an rf discharge [0. The discharge 
is obtained by applying pulsed rf power to a helical res¬ 
onator (HNMR) located close to the QM and operating 
at Ri 910 MHz. An extra chamber with another RF coil 
attached to the sample cell from the top was used as 
atomic hydrogen gas source. Atoms in the gas phase 
can be stabilized if the sample cell walls are covered by 
superfluid helium film. Having a possibility of filling the 
sample cell with the gas of H or D provides good magnetic 
field markers for the accurate measurement of the shifts 
of the ESR lines originating from the atoms in solid hy¬ 
drogen films. Measuring heat released in recombination 
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Figure 2. Dependence of H and D atom concentrations in 
solid HD on time during the process of atom accumulation 
by running discharge in the sample cell and after discharge 
was switched off. The D concentration reaches steady state 
in a time comparable with the time constant of reaction 
2 X 10“^ s while H concentration continues to grow even further. 
The D concentration rapidly decreases due to reaction @ 
after the discharge is switched off. 



Figure 3. Deuterium atom concentration decay in solid HD 
at different temperatures. Note that the slopes of the lines 
are almost parallel. 


of gas phase atoms allows independent calibration of the 
absolute number of spins in a sample. For preparation of 
our samples we used hydrogen isotopes purchased from 
Linde AG, D 2 gas with 0.23% residual concentration of 
HD, and HD gas with 3% of H 2 . In order to obtain the 
matrix molecules in their ground rotational state, prior 
to the condensation of the D 2 gas into the sampk cell we 
kept it overnight in the para-ortho converter [l7|. This 
ensures high concentration,~90%, of ortho-deuterium in 
our samples. 

III. EXPERIMENTAL RESULTS 

The experimental study of the exchange reaction (|3]) 
was carried out in two different matrices: a pure HD ma¬ 
trix and a D 2 matrix with a 0.23% HD admixture. These 
two systems allowed us to study two different approaches: 
the first case, where the D atoms are surrounded by 12 
HD molecules, no D-atom diffusion is needed for the re¬ 
action to proceed. In the latter case, D atoms need to 
diffuse a distance of several lattice constants before they 
will encounter a HD molecule and accomplish the ex¬ 
change reaction. 

At the beginning of the measurements with the pure 
HD sample, a 200 nm HD film was deposited onto the 
quartz microbalance. Then the rf discharge was started 
in order to begin accumulating H and D atoms in the 
HD solid. The sample cell temperature was increased 
to 0.7 K due to the discharge. During accumulation, H 
and D atom populations were monitored by ESR (Fig. 
@) and when the D concentration saturated, the dis¬ 
charge was turned off. Typically the discharge was run 
for 5-12 h. The maximum H and D concentrations we 
achieved were 2 x 10^®cm“^ and 1.5 x 10“^®cm“^, re¬ 


spectively. Then the system was set to a desired tem¬ 
perature and the measurements of the D atom decay be¬ 
gan. ESR spectrum of atomic D contains three lines sep¬ 
arated by 78 G due to hyperfine interaction[l^ . In ther¬ 
mal equilibrium the lines have equal area, and therefore 
for monitoring the D concentration evolution is sufficient 
to record only one line. This has been done every ten 
minutes until the D lines vanished in the noise. After 
finishing one measurement the process of atom accumu¬ 
lation was repeated. Even after several accumulations 
we did not observe differences in the sample properties. 
Moreover exchange reaction rates measured for different 
HD films were quite reproducible. The D atom decay 
curves obtained at different temperatures are presented 
in FigO The rather high temperature during discharge 
(0.7 K) provided a problem for the lowest temperature 
measurements because the system required 1-2 h to cool 
down to T« 130 mK during which the D concentration 
had already decreased considerably. 

Absolute concentrations of atoms were determined via 
the effect of broadening the ESR lines due to dipole- 
dipole interaction between the atoms[l3|. At the high 
densities considered in this work this broadening mech¬ 
anism is substantially stronger than other effects related 
with the nuclear dipolar moments of the matrix, and 
provides reliable calibration of the ESR line area ver¬ 
sus concentration of atoms. Independent verification of 
this calibration was performed calorimetrically, measur¬ 
ing the heat released in recombination of H atoms in the 
gas phase. We evaluate possible systematic error in the 
concentration measurement as 20%. 

The reaction ([4]) is described by the kinetic equation 

= -fee. [HD] [D], (5) 

where kex is the rate constant. Because [D] <C [HD], [HD] 






4 



i 

EH-HD^ D 2 +H 


D+D-^'Eh 

0 

m 

0 

[3] 


pure HD 


D2:0.23%HD 

□ 

[8] 

A 

[13] 

5 

EH-HD^ Eh+H 



n 

[5] 

< 

[24] 


E)2:0.23%HD 

0 

[20] 

0 

[14] 

V 

[25] 


Figure 4. The H in H 2 recombination rates (green symbols), D 
in D 2 recombination rates (red symbols) and the rates of the 
exchange reaction ([4]) (blue symbols) as a function of temper¬ 
ature. The rates of the exchange reaction ((4| in D2:0.23%HD 
matrix are shown as yellow triangles. The filled symbols cor¬ 
respond to experimental data obtained in this work. The lines 
are shown only as a guide. 

can be considered constant and eq. (O becomes a pseudo 
first-order reaction. Then slopes of the linear fits of semi¬ 
log plot in Fig. [3] give the value of k [HD] where we used 
the density of HD molecules equal to 0.0486 moP^cm“^ 
[13. This measurement is free of errors in the absolute 
concentration calibration because the first-order reaction 
rate depends only on the time constant. Possible recom¬ 
bination of atoms is negligibly small since the change of 
the total concentration [H] -|- [D] was not observed within 
the accuracy of the measurements. Our results for kex 
and values of these rate constants from the literature are 
presented in FigU) Takwanagi et al. carried out theoret¬ 
ical calculations of fcea;[^. Prior to this work the lowest 
temperature measurement was done by Kiselev et al. (l9| 
and Bernard et al. who studied the reaction (|T1) in 
HD-D 2 nanoclusters of impurity helium condensates at 
temperature of 1.35 K. The rates obtained in these works 
are in a good agreement with our results. 

One may expect that electrons and ions resulting from 
the discharge may also be trapped in our solid samples. 
Sometimes we observed an extra line in the center of 
ESR spectrum attributed to the trapped electrons (2l| . 
Presence of ions was also observed in the studies of HD 
samples irradiated by 7 —and cc-rays by Kumada et 
al. Although such impurities may distort the lattice and 
severely influence its properties, we do not think that this 
has any importance in the measurements of the exchange 
reaction in pure HD, since no any diffusion of reagents is 
required for the reaction. We have not found differences 
in the H or D atom properties in samples which exhibit 
an ESR line at the free electron resonance pH . 



Tigure 5. The evolution of H and D atom fractions in a D 2 
matrix with a 0.23% HD admixture as a result of reaction 
(|4|l at temperature 150 mK (grey diamonds) and 1.35 K (blue 
circles). The total concentration of H and D atoms measured 
at 150mK is shown by red squares. The [H],[D]/([H] + [D]) 
axis is shown broken to emphasize the growth of the H atom 
fraction. 

In a second part of our work we studied the exchange 
reactions in D 2 sample prepared from the pure D 2 gas, 
the most pure commercially available: with an HD ad¬ 
mixture of 0.23% (6.9 X 10^®cm“^) as specified by the 
manufacturer (Linde AG). It turned out that due to the 
exchange reaction ([4]) even this tiny amount of HD re¬ 
sults in a rather large concentrations of H atoms in the 
matrix. 

A 200 nm D 2 film was condensed directly from the 
room temperature reservoir similar to that described be¬ 
fore. The measurement was started after accumulating 
atoms in the film by running a discharge in the sample 
cell for two days and subsequently cooling to T=150mK 
(FiglS]). The concentrations of D and H atoms ap¬ 
proached 2.3xl0^®cm“^ and 2 . 1 xl 0 ^®cm“^ after ^ 200 h 
of measurement. The H:D ratio is much greater than 
what one should expect from direct dissociation of H and 
D atoms during discharge and assuming that no exchange 
reaction was taking place. Based on the gas composition 
the D:H ratio should be about 1700:1. The total H and 
D atom concentration remained constant in the course of 
measurement when discharge was kept off which means 
that neither H nor D recombination was taking place. 
The concentrations of the reactants D and HD in the 
D2+0.23% HD sample were of the same order during the 
measurement. Therefore, for determination of the reac¬ 
tion rates we need to solve the equation ([5]) assuming 
true second order process. The solution can be presented 
in a form 

1 , JHD][D]o _ 

[HD]o-[D]o [D][HD]o “ 

where [HD]o and [D]o are the concentrations at the begin¬ 
ning of measurement while [HD] and [D] are the time de- 
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pendent concentrations of the reactants during the course 
of measurement. Concentration of D atoms we deter¬ 
mined from the ESR lines area using the absolute density 
calibration as described above. The values of [HD] were 
estimated indirectly from the [H] growth as [HD]=[HD]o- 
[H], assuming that the exchange reaction is the only pro¬ 
cess leading to the decrease of [HD]. Using experimental 
data in Figl5| we determined the effective exchange re¬ 
action rate at 150 mK kexD 2 = 9(4) x 10“^® cm®s“^. 
We repeated the described above measurement of the ex¬ 
change reaction at temperature of 1.35 K and observed 
practically same time dependence, as can be seen in Fig. 
®. Unfortunately, the long-term stability of our ESR 
spectrometer degraded at temperatures above 1 K, and 
was not sufficient for accurate monitoring of the H and D 
densities. Therefore, in FiglSjwe plotted the fractional ra¬ 
tios of densities [H]/([H]+[D[) and [D[ / ([H[-|-[D[), which 
were not sensitive to the slow drifts of the spectrometer 
sensitivity and better represent the changes caused by the 
exchange reactions. Using the 1.35 K data in Fig. jSjwe 
extracted the rate constant kexD 2 ~ 8 X 10 cm®s ^ 
very close to the value obtained at 150 mK. Such be¬ 
haviour indicates that the exchange reaction proceeded 
with the same rate having very weak or no tempera¬ 
ture dependence. We attempted similar measurements 
at higher temperatures of 2.5 and 3.5 K. In addition to 
the conversion of D into H we also observed a decrease of 
the total density of atoms caused by their recombination. 
Since the recombination includes several channels includ¬ 
ing H-l-H, H-|-D, and D-|-D processes, and also because of 
the increased drifts at higher temperatures, we have not 
attempted to extract any quantitative data from these 
measurements. Such work requires optimization of the 
experimental setup for operation at higher temperatures 
and may be performed in the future. 

IV. DISCUSSION 

In this work we performed measurements of the hydro¬ 
gen isotope exchange reactions in the matrices of pure HD 
and D 2 . Other well known reactions which occur in these 
systems are the reactions of recombination of atoms into 
molecules. These reactions are also of the second order in 
density. This allows direct comparison of the rates, which 
we perform in Fig |4| where the rate constants of the ex¬ 
change reaction (|4]), the recombination rate constant of 
H in H 2 , and D in D 2 are presented for the temperature 
range 0.15-4 K. The plot includes results of the present 
work as well as a compilation of our previous results and 
results of other groups. The scarcity of the reactants ([H] 
in H 2 [HD[ in pure HD) requires a much longer time to 
detect the H atom decay although the reaction rate turns 
out to be greater than that of the exchange reaction (jj]) 
(FiglU). 

The measurements in pure HD demonstrate that the 


exchange reaction takes place at the lowest temperature 
at the same rate as at 1.5 K without any clear temper¬ 
ature dependence. Such behaviour is expected for tun¬ 
nelling reaction and means that the chemical reactions 
of this type may proceed down to absolute zero of tem¬ 
perature. In this case we studied the pure exchange re¬ 
action o, because each D atom is surrounded by 12 HD 
molecules needed for the reaction, and diffusion stage is 
not required. The total H and D atom concentration 
in the HD samples we studied was about 3xI0^®cm“® 
which corresponds to the mean distance of 3 nm or 8 
lattice constants between them. The H-D or D-D in¬ 
teraction at such distances becomes negligibly small and 
cannot influence the rate of the exchange reaction (jj]) 

Q. 

In the samples of D 2 with small concentration of HD 
dissociation of atoms by the discharge produces mainly 
D atoms which on the average are located quite far away 
from HD impurities. We may neglect production of H 
and D atoms by the direct dissociation of HD, because 
its concentration in the matrix is much smaller than D 2 . 
Since each D atom has 12 neighbours, the probability 
that one of the neighbours is HD may be estimated as 
1 — 0.9977^^ ^ 2.7%. However, as we can see from Fig. 
[51 almost 10% of D atoms were converted to H which 
provides evidence for D atom diffusion. The complete 
sequence of the D-to-H conversion in this case includes 
the diffusion process D -|- D 2 D 2 -I- D and the exchange 
reaction itself. This provides an opposite limit to 
the described above experiments with pure HD where 
diffusion is not required at all. 

Absence of temperature dependence in our data on the 
exchange reaction rate in the D 2 matrix indicate that the 
process is not diffusion limited, and the reaction itself is 
the limiting stage. Then, comparing the absolute values 
of the rate constant we note that in the HD matrix it 
is only a factor of 3 larger than in D 2 . Since the diffu¬ 
sion stage could further slow down the effective rate of 
the exchange reaction in the D 2 matrix, the difference in 
the tunnelling rates of the reaction may be even smaller. 
Considering that in hep lattice of pure HD each D impu¬ 
rity is surrounded by 12 neighbours, one may expect that 
in D 2 matrix the exchange reaction should be at least 12 
times slower. We would like to point out also on a sub¬ 
stantial difference in the initial and final states of the 
reactants in these two matrices. For the exchange reac¬ 
tion o in pure HD the initial state has a single crystal 
defect, which can be represented in I dimensional case 
as ...HD-HD-D-HD-HD... After the act of the exchange 
reaction we obtain two defects: ...HD-HD-D 2 -H-HD..., 
which in total create much larger distortion of the lat¬ 
tice than in the initial state. On the contrary, in the 
D 2 matrix we have initially ...D 2 -D 2 -D-HD-D 2 ..., which 
after the exchange is converted to: ...D 2 -D 2 -H-D 2 -D 2 ... 
Now the situation is opposite, two lattice defects are con¬ 
verted into one. Although, the mass is conserved in the 
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exchange reactions, the volume of the lattice distortion 
is much larger in the case if two neighbouring defects, 
which may lead to a larger difference in the energies of 
the initial and final states for the reaction in the HD ma¬ 
trix. We consider this as a possible explanation for the 
observed difference in the exchange reaction probabilities 
in HD and D 2 matrices. 

The recombination of H atoms in H 2 solids also in¬ 
cludes two stages: approaching of H atoms towards each 
other by a distance of one lattice constant by a series of 
tunnelling reactions H-I-H 2 ^ H 2 -I-H followed by the for¬ 
mation of a H 2 molecule by the recombination reaction 
H -|- H —>■ H 2 with a transfer of recombination energy to 
the lattice. For spin-polarized atoms in the lattice there 
appears to be several fast depolarization processes which 
lead to a flip of the electron spin of one of the atoms and 
subsequent recombination. It has been shown that the 
rate of this process does not depend on magnetic field at 
temperatures above 1.3 K [^, and the diffusion stage is 
always the limiting stage of the recombination. 

Hydrogen migration in a perfect H 2 crystal proceeds 
via resonant tunnelling when the initial and final states 
of a tunnelling event coincide. The presence of H atoms 
(or other crystal imperfections) in the lattice causes the 
mismatch of the initial and final levels which greatly re¬ 
duces the tunnelling rate. The mismatch increases when 
H atoms approach each other and finally it prevents them 
from occupying the neighbouring lattice sites where they 
are able to recombine. The level mismatch can be over¬ 
come by phonon assistance via a single-phonon direct 
process or a two-phonon Raman process, which lead to 
a strong temperature dependence of the recombination 
rate. 

Data for the recombination of H atoms in solid H 2 
summarized in Fig. |4] showed a strong temperature de¬ 
pendence with two distinct ranges: 1.3-4 K where it de¬ 
pends linearly on temperature and below 1K where the 
rate constant decreases substantially and become immea¬ 
surably small at 150 mK[l^. Recombination becomes 
strongly inhibited at the lowest temperatures, about 
150 mK where no phonons are available to compensate 
for the level mismatch. Ahokas et al. showed that H- 
atom recombination can be initiated at temperatures be¬ 
low 500 mK if phonons are injected into the system [l^ 
as a result of recombination of the gas-phase Hj, atoms 
on the H 2 film surface. 

Recombination of atomic deuterium in solid D 2 has 
similar dependence on temperature, and is also diffusion 
limited. The absolute values of the rate constants are 
two orders of magnitude smaller, and approaching to 1 
K one cannot detect any recombination making its mea¬ 
surement impossible. The upper limit estimate for the 
D atom recombination constant at 150 mK found from 
Figl5]is equal to krD 2 — 3(2) x 10“^®cm“^. The value of 
krD 2 T=150mK is much smaller than that measured 
at 4.2 K and 1.5 K reported earlier (Fig]?]) 241, [IHl which 


gives evidence for a strong temperature dependence of 
the recombination constant of D atoms below 1.5 K, sim¬ 
ilar to that observed for H in H 2 below 1 K [3| ■ Heavier D 
atoms in a D 2 crystal have smaller tunneling probability 
than for H in H 2 and diffuse much slower. 

As we see, the diffusion rate of the reactants of the 
exchange reaction (|3]) in the D 2 matrix is substantially 
faster than the diffusion of two D atoms approaching each 
other for recombination. We suggest two possible reasons 
for that. First, there is a smaller mass difference Am be¬ 
tween a HD moleeule and hosting D 2 molecules relative 
to the mass M of the host molecule: Am/M = 1/4 com¬ 
pared to that of a Am/M = 1/2 for D in D 2 (and also 
for H in H 2 ). This implies a smaller perturbation of the 
periodic potential of the matrix in the former case. Sec¬ 
ond, it is known that the HD molecules in their ground 
state may rapidly move in the D 2 crystals via tunnelling 
exchange with ortho-D2[2^. This process may go even 
faster than the diffusion of D atoms due to the reaction 
(p|). In this case, the energy level mismatch should be 
overcome as well. As a result, it becomes easier to ap¬ 
proach each other for D atom and HD molecule than for 
two D atoms. The recombination rate of D atoms in a 
D 2 matrix should be slower than the rate of the exchange 
reaction (|4]) as observed in the present work. 

The exchange reaction ([3]) is expected to proceed two 
orders of magnitude faster than the reaction (|4]) and can¬ 
not be measured in a pure H 2 matrix with the accumu¬ 
lation technique used in our experiments. Measurements 
of samples with a small fraction of H 2 are diffusion lim¬ 
ited and therefore do not represent the true reaction rate 
[3 . But such measurements will yield information on the 
diffusion of D if the reaction ([3]) remains fast at ultralow 
temperatures. One should be able to make lower limit 
estimates of the reaction ([3]) by comparing the yields of 
D and H in samples with different D 2 : H 2 mixtures. For 
H 2 fractions above 25% no D lines are expected [ 2 ^. 


V. CONCLUSIONS 

In conclusion, we have reported on the first experi¬ 
mental observation of the isotopic exchange tunnelling 
reaction, D -|- HD ^ D 2 -I- H, taking place in a solid ma¬ 
trix of hydrogen deuteride at temperatures below 1K. 
The reaction rate, kex, was measured within the tem¬ 
perature range 1.5 K - 130 mK where it was found to be 
nearly temperature independent, k^x ~ 3 x 10“^^cm^s“^. 
The exchange reaction also takes place at T=150 mK and 
1.35 K in a D 2 matrix with 0.23% HD admixture where D 
atoms need to migrate by the distance of several lattice 
constants in order to participate in the reaction. The re¬ 
action rate did not change upon raising temperature from 
150 mK to 1.35 K. The fact that the isotopic exchange re¬ 
action D-|-HD^D 2 +H proceeds there suggests that the 
exchange reaction D-I-D 2 —>'D 2 -t-D which governs D atom 
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diffusion in a D 2 matrix is substantially faster than the 
exchange reaction D + HD —D 2 + H. Reconibiiration of 
D atoms becomes suppressed due to the energy level mis¬ 
match appearing when two atoms approach each other. 
The level mismatch depends on the mass difference be¬ 
tween the atomic impurity and the molecule of the host 
lattice, and is substantially smaller for the case of HD in 
the D 2 matrix. 

These observations have a bearing on the rate of the 
exchange reaction H-I-H 2 ^H 2 +H also does not depend 
on temperature if H and H 2 occupy ireighboring sites 
even though the rate of the atomic hydrogen recombina¬ 
tion reaction, H -|- H —H 2 , is substantially reduced upon 
lowering temperature below 1 K and it almost completely 
vanishes at 150 mK. Therefore we anticipate that recom¬ 
bination of H atoms in solid H 2 is mainly limited by the 
diffusion stage, the exchange reaction H -f H 2 —H 2 -t- H, 
which is inhibited in samples with high atomic concentra¬ 
tion at low temperatures where the condition of resonant 
tunneling becomes violated. 

The experiments show strong decreasing of recombi¬ 
nation rates of H (D) atoms in solid H 2 (D 2 ) due to 
suppression of diffusion of atoms at temperatures below 
1 K. New theoretical approaches are needed to under¬ 
stand this phenomenon. 
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